Directed anomalous diffusion without a biased field: A ratchet accelerator 



^1" 
O 

o 

(N 
>v 



(N 

Q 

U 
d 



m 
o 
in 
o 

o 



> 
X 



Jiangbin Gong Q and Paul Brumer 
Chemical Physics Theory Group, Department of Chemistry, 
University of Toronto, Toronto, Canada M5S 3H6 
(Dated: February 8, 2008) 

Directed classical current that increases linearly with time without using a biased external field 
is obtained in a simple model Hamiltonian system derived from a modified kicked rotor model, by 
breaking the spatial symmetry of the transporting regular islands in the classical phase space. A 
parallel study of the corresponding quantum dynamics suggests that although quantum coherence 
effects suppress the directed current and can even induce a current reversal, directed quantum 
current that increases linearly with time can nonetheless be realized in a quantum system that is 
far from the classical limit. 
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I. INTRODUCTION 

Ratchet effects, i.e., directed transport without a bi- 
ased field due to a broken spatial-temporal symmetry, 
are of great interest to the understanding of molecular 
motors in biological systems and have been under in- 
tense investigation in a variety of simple model systems 
0, 0] . Early studies focused on the role of external noise, 
which was then replaced bydeterministic chaotic dynam- 
ics with dissipation 0,0, IE 01- Recently, it was shown 
that even purely Hamiltonian dynamics with both reg- 
ular and chaotic phase space structures [3, la, |3, llfj or 
with complete chaos 0, 0, 0, is capable of gener- 
ating ratchet effects. As such, there has been a great 
interest in studying simple model Hamiltonian systems 
in an effort to observe and understand ratchet effects. 

The kicked rotor (KR) model 01 and its variants 
are ideal systems for such studies. For example, delta- 
kicked systems with an alternating kicking period and an 
asymmetric potential [9|, llll Il2j | , or with multiple exter- 
nal fields [l(| have been considered in connection with 
ratchet effects. These studies are stimulating since the 
classical dynamics of the KR (i.e., the standard map) 
and the associated quantum dynamics have long served 
as paradigms for classical and quantum chaos |lj|- In- 
deed, the KR has been realized in several atom optics 
experimental groups 0] , and is also of considerable in- 
terest in other fields such as condensed matter physics 
HHm, molecular physics pl.Hll20ll2l|. and quantum 
information science [23, [23 ■ 

Although not often thought of as such, the basic idea 
that underlies the above-mentioned studies of ratchet ef- 
fects in the KR, i.e., to manipulate symmetry properties 
to induce a directed current, is not new from a coherent 
control perspective. Indeed, fifteen years ago Kurizki, 
Shapiro, and Brumer |24j proposed a simple coherent 
control scenario that can generate photocurrents in semi- 
conductors without a biased voltage. This proposal and 
its various extensions [25l |2(j have been experimentally 
realized [H |H |29| . 

In this paper we introduce a totally different type of 
Hamiltonian ratchet effect in a simple delta-kicked model 



system, which we call a "ratchet accelerator". Without 
a net external force (i.e., the external force averaged over 
space is zero), the ratchet accelerator can generate di- 
rected classical anomalous transport, with the net cur- 
rent (defined below) that accelerates linearly with time. 
This is made possible by taking advantage of the large 
transporting islands in the classical phase space of a 
modified kicked rotor system recently proposed by Gong, 
Worner, and Brumer |30ll31| . in conjunction with a sec- 
ond kicking field that breaks the spatial symmetry. The 
corresponding quantum dynamics is also studied, shed- 
ding considerable light on quantum-classical differences 
and correspondence in the presence of anomalous diffu- 
sion. For example, we show that due to the tunneling 
between a chaotic sea and a relatively large regular re- 
gion of the classical phase space, the quantum directed 
current can be in the opposite direction of the classical 
directed current. 

This paper is organized as follows. In Sec. II we de- 
scribe the background and the motivation of our model 
as a ratchet accelerator. The classical results of directed 
anomalous diffusion are presented in Sec. Ill, followed 
by the analogous quantum results in Sec. IV. Conclud- 
ing remarks are given in Sec. V. 



II. THE MODEL SYSTEM 

Consider the following model Hamiltonian 
H(p,x,t) = p 2 /2m + AiCOs(:e/:eo)^ f(n)S(t/T - n) 

n 

-\ 2 sm(2x/x )J2 6 ( t / T - n )> C 1 ) 

n 

where p is the momentum, m is the particle mass, x is the 
conjugate position, Ai and A2 are the strength of the first 
and second kicking fields, T is the time interval between 
kicks, and the function f(n) will be determined below. 
The associated classical map is given by 

Pn = pat-i + Ki/(A r )sin(x JV -i) + 2k 2 cos(2£at_i); 
x N = x N -i + Pn, (2) 
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where x = x/x° is the scaled dimensionless position vari- 
able, p = pT '/ '(mio) is the scaled dimensionless momen- 
tum variable, K\ = \\T 2 / (mip) , K2 = A2T 2 / '(mxg) , and 
(Pn,xn) represents the phase space location of a clas- 
sical trajectory at (N + 1 — + )T. The corresponding- 
quantum unitary evolution operator for propagating from 
time (N - 0+)T to time (N + l- 0+)T is given by 
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exp[—if(N)ki cos(x)] exp[ifc 2 sin(2x)],(3) 



with dimensionless parameters k% = XiT/h, fc 2 = X2T /%, 
t = hT/mx 2 ,. Let p = pT/(mxo), where p — —ihd/dx, 
then one has [x,p] = it, and p — —ird/dx, indicating 
that r plays the role of an effective Planck constant. Note 
that K\ = k\T and «2 = k%T. 

Clearly, if A 2 = and f(n) = 1 then Eq. Q reduces 
to the standard kicked rotor model 0] and the classical 
map of Eq. becomes the standard map. If A2 = 0, 
f(n) = g{n), where g(n) = 1 for n = Aj + 1, Aj + 2, and 
g(n) = — 1 for n = 4j + 3, Aj + A, where j an inte ger, then 
the system becomes that in our previous work |30l l3l|. 
Our accelerator, introduced here, assumes A2 7^ and 
f(n) — g(n). That is, the system is kicked by two fields, 
one of which reverses the sign of its kicking potential after 
every two kicks. 

Note first the existence of "accelerating trajectories" 
in the classical delta-kicked dynamics. In particular, 
for particular values of Hi, the standard map [K2 = 0, 
f(n) = 1] can generate trajectories whose momentum in- 
creases (or decreases) linearly with time (at least on the 
average). These trajectories are a class of transporting 
trajectories To see this consider the initial condi- 
tions: (p — 2ttIx,x — ±7r/2) for k = 27rZ 2 , where l\ 
and I2 are integers. Clearly, these phase space points 
are shifted by a constant value (±27^2) in p after each 
iteration, resulting in a quadratic increase of rotational 
energy. These accelerating trajectories are rather stable 
insofar as they may persist for values of k close to 27r/ 2 
(with their average momentum shift after each iteration 
oscillating around the constant value ±271^2), thus giv- 
ing rise to transporting regular islands ||| , also called the 
"accelerator modes" in the standard map case j33, ■ 
Dramatically, if f(n) = g{n) and K2 = 0, then there 
can exist much larger transporting islands [3(j. These 
islands are associated with the marginally stable points 
(p = (2Zi + 1)tt,x = ±7r/2) for k = (2l 2 + 1)tt. Trajecto- 
ries on these islands will be shifted by a constant value 
[±(2?2 + 1)tt] in p after each kick. 

If classical trajectories are launched from the trans- 
porting regular islands, they simply jump to other similar 
islands located in adjacent phase space cells. For trajec- 
tories initially outside the transporting regular islands, 
the "stickiness" of the boundary between the islands 
and the chaotic sea induces anomalous diffusion over the 
phase space. For example, the square of the variance 
of momentum increases nonlinearly, but not quadrati- 
cally. This is intrinsically different from the case of nor- 



mal chaotic diffusion in which the square of the variance 
of momentum increases linearly with the number of kicks. 

As in previous studies on delta-kicked systems, one 
needs to break the spatial symmetry of the kicking po- 
tential to induce ratchet effects. More specifically, our 
choice of the above model Hamiltonian is motivated by 
the desire to have large transporting islands in the clas- 
sical phase space |30| and to break the spatial symme- 
try of the transporting regular islands. As shown below, 
the symmetry can be indeed broken by a second kicking 
field, and yet the remaining transporting islands are still 
of a significant size and have unidirectional transport- 
ing properties. This being the case, the mechanism of 
the Hamiltonian ratchet effects described below becomes 
transparent without using previous analyses of the re- 
lationship between broken spatial-temporal symmetries 
and ratchet effects H IH HI IH . 



III. CLASSICAL RESULTS 

To demonstrate the feasibility of our basic idea, in 
this section we present our results using some specific 
computational examples. Figure^ displays the classical 
phase space structure of the standard map with Ki = 3.8, 
K2 = 0, and f(n) — 1. The regular islands in Fig. are 
not transporting since the momentum of the trajecto- 
ries launched from these regular islands is bounded and 
oscillates periodically. Figure^) shows that if the first 
kicking field reverses its potential after every two kicks, 
then the regular non-transporting islands in the stan- 
dard map are destroyed and new regular islands emerge. 
A simple computation reveals that those islands seen in 
Fig. are transporting. In particular, classical tra- 
jectories launched from regular islands on the left side 
(x < 7r) of the phase space cell will accelerate, with their 
p increased by ~ 7r after each kick. Due to the spa- 
tial symmetry, these islands always have partners with 
the opposite transporting property: trajectories launched 
from those islands on the right side (x > ir) of the same 
phase space cell will accelerate in the opposite direction, 
with their p decreased by about 7r after each kick. These 
properties of the transporting regular islands seen in Fig. 
^> indicate that they are associated with the marginal 
stable points p = it, x = ±7r/2 in the case of K\ = tt. 
Dramatically, upon introducing the second kicking field 
in Eq. 0, e.g., A2/A1 = k 2 /ki = 0.2/3.8, the trans- 
porting islands no longer appear in pairs as in Fig. ^3. 
Rather, as seen in Fig. the spatial symmetry of the 
phase space structures is clearly broken, but the size of 
the main transporting islands that persist remains signif- 
icant. Specifically, all the trajectories launched from the 
regular structures clearly seen in Fig. ^ wm still accel- 
erate linearly with time; the remaining part of the phase 
space is exclusively occupied by the chaotic sea. That 
is, there does not exist similar transporting islands that 
decreases p in a linear fashion. 

For a classical ensemble that uniformly covers the en- 
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kicking fields is zero and the initial classical ensemble is 
both spatially symmetric and time-reversal invariant, all 
the trajectories tend to be accelerated in one direction. 
Note that, by changing the strength of the second kick- 
ing field or varying the relative phase between the two 
kicking fields (e.g., the second kicking field is replaced by 
A2sin(2a;/xo + P)^2 n ${t/T — n), where (3 is a relative 
phase factor), one can also obtain a r = —tt. Hence, the 
acceleration can be also obtained in the reverse direction. 

For an initial classical ensemble with p = and x uni- 
formly distributed from — oo to +oo, Fig. [2b displays the 
time dependence of the classical net current (p), where (•) 
denotes the ensemble average. It is indeed seen that the 
classical net current is in the negative direction, and ac- 
celerates linearly with time. The acceleration rate is also 
consistent with Eq. (0}. Figure [2b displays a log — log 
plot of the time dependence of the associated momen- 

1/2 

turn variance Ap = ((p 2 ) — (p) 2 ) . A linear fit of the 
result in Fig. [2b gives Ap ~ ]\[0-65^ wnere ]\[ [ s the num- 
ber of kicks. This confirms that the diffusion dynamics 
is anomalous and immediately gives Ap/(p) ~ 1/JV - 35 . 
Hence, as far as the first and second order statistical mo- 
ments are concerned, the momentum fluctuations of the 
classical ensemble become less and less important as time 
increases. In this sense, for sufficiently large times, al- 
most all trajectories in the classical ensemble are accel- 
erated in the negative direction. 



IV. QUANTUM RESULTS 



FIG. 1: Classical phase space structures of (a) the standard 
map [«i = 3.8, f(n) = 1, Ki = 0], (b) a modified kicked rotor 
model [ki = 3.8, f{n) — g(n), K2 = 0], and (c) a ratchet 
accelerator [m — 3.8, f(n) = g(n), k,2 = 0.2]. All variables 
are in dimensionless units. Note that the regular islands seen 
in (b) and (c) are transporting while those in (a) are not. 

tire phase space cell shown in Fig. ^j, the spatial and 
temporal periodicity of the system requires that the av- 
erage acceleration rate is zero [lH, i.e., 

a c A c + a r A r = 0, (4) 

where A c and A r denote the areas of the chaotic sea 
and the regular transporting region, respectively, with 
their associated average acceleration rate denoted by a c 
and a r . Consider now an initial classical ensemble with 
p = and i uniformly distributed between — oo and +oo. 
Such an ensemble is easy to simulate both classically and 
quantum mechanically. Since this initial ensemble is en- 
tirely located in the chaotic sea (see Fig. QJ;), Eq. (0J 
predicts that the average acceleration rate of this classi- 
cal ensemble is given by — a r A r /A c . This also indicates 
that the larger the transporting island is, the faster the 
directed acceleration would be. In the case of Fig. QJ;, 
a r = it, so more trajectories travel to the regime of neg- 
ative x. That is, even though the average force of the 



How then does this ratchet accelerator manifests it- 
self in the corresponding quantum system and how small 
the effective Planck constant should be in order to ob- 
serve the analogous ratchet effects in the quantum dy- 
namics? To carry out the quantum calculations we im- 
pose a periodic boundary condition on the system, i.e., 
\ip(x)) = \ip{x + 2tt)), where \iJj) is the wavefunction of 
the system. This periodic boundary condition suggests 
that St is now to be understood as an angle variable, and 
is consistent with the fact that the classical dynamics is 
invariant if x — ► x + 2tt. Further, our previous work |30| 
suggests that this choice of boundary conditions should 
not affect the essence of the quantum dynamics. 

Figure [3 shows the net quantum current, i.e., the ex- 
pectation value of p for two values of the effective Planck 
constant r, with the initial state given by \ip(t = 0)) = 
|0). Comparing the results here with that in Fig. [2b, one 
sees that the magnitude of the quantum directed cur- 
rent can be much smaller than that of the classical di- 
rected current. Interestingly, the directed current in the 
quantum case displays much more complicated behavior 
than the linear acceleration seen classically (some slow 
oscillatory behavior of the quantum current over much 
longer time scales is also observed, but should be of much 
less experimental interest). This can be qualitatively ex- 
plained in terms of a combination of a dynamical local- 
ization effect |l4j and the quantum tunneling between 
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FIG. 3: The time dependence of the quantum dimensionless 
directed current defined as the expectation value of p for two 
values of the effective Planck constant r, with the initial con- 
dition given by |0) and other parameters are the same as in 
Fig. Results here should be compared with those in Fig. 

m 



This means that for experiments using cold sodium atoms 
|l5j | the directed and relatively large quantum current ob- 
tained with 2000 kicks would correspond to an average 
drift velocity on the order of 1 m/s. Note however, this 
value of t is about one order of magnitude smaller than 
what is being examined in atom optics experiments |15| . 
but may be achievable in the future. 



FIG. 2: (a) The linear time dependence of the classical di- 
rected current defined as the ensemble averaged value of the 
momentum, with the initial condition p = and x uniformly 
distributed between — oo and +oo. (b) The associated time 
dependence of the momentum variance in a log, log plot. The 
system parameters are the same as those in the case of Fig. 
All variables are in dimensionless units. 



the chaotic sea and the transporting islands 30, 37]] ■ In 
particular, since the dynamical localization in quantum 
delta-kicked systems always suppresses the diffusion, the 
quantum directed current is not expected to grow all the 
time. More importantly, since the chaotic sea and the 
transporting islands have opposite acceleration signs [see 
Eq. , the quantum tunneling between them |30|] nec- 
essarily decreases the directed current. Indeed, in the 
case of t = 0.2, one sees that for most of the time the 
quantum current is in the opposite direction of the clas- 
sical current. Note that, unlike a previous observation of 
tunneling-induced current reversal in quantum Brownian 
motion [38|, here the tunneling is between a chaotic sea 
and a regular region and the system is totally isolated 
from the environment. 

The results in Fig. |21 also suggest that t = 0.1 is 
small enough to observe the signature of a ratchet accel- 
erator in a quantum delta-kicked system. For example, 
the quantum current for r = 0.1 does increase almost lin- 
early for as many as 2000 kicks, with the average increase 
rate about three times smaller than in the classical case. 




-150 -100 -50 50 100 150 
m/1000 



FIG. 4: The probability P(m) of finding the quantum system 
in state \m) after 4000 kicks, with the initial state given by 
|0). Solid line is for r = 0.1 and dashed line is for r = 0.2. 
Other parameters are the same as those in the case of Fig. 



It is also interesting to examine the lineshape of dy- 
namical localization with the spatial symmetry of the 
transporting regular islands broken by the second kick- 
ing field. Figure 0] displays the associated distribution 
function P(m) after 4000 kicks for the same two values 
of t. It is clearly seen that for both cases P(m) is asym- 
metric. Due to the stickiness of the boundary between 
the transporting islands and the chaotic sea and the as- 
sociated quantum tunneling, the probability of finding 
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the system in state |m) with a large and positive m can 
be much higher than that of finding the system in state 
— m). This does not contradict the previous observa- 
tion that after an initial stage the quantum current for 
t = 0.1 is negative, because the tail of the statistical 
distribution function P(m) has little effect on the net 
current, the statistical moment of the lowest order. The 
huge difference in P(m) between the two cases shown in 
Fig. Q]with slightly different values of r is also notewor- 
thy. This extreme sensitivity of P(m) to the value of r 
is more or less related to the fact that the values of r 
we consider are quite comparable to the area of the main 
transporting island seen in Fig. 

V. CONCLUDING REMARKS 

Using a simple model Hamiltonian system we have 
demonstrated a new ratchet effect by considering directed 
classical anomalous diffusion. We have shown that it is 
possible to achieve directed acceleration with a fixed ac- 
celeration rate with two kicking fields that have zero av- 
erage gradient. This is the case even though the initial 
classical ensemble is both spatially symmetric and time- 
reversal invariant. We have also presented some compu- 
tational examples of the quantum dynamics with differ- 
ent values of the effective Planck constant, ft is found 
that quantum coherence effects may induce a reversal of 
directed current. It is also found that directed quantum 
current that accelerates linearly with time without us- 
ing a biased external field should also be achievable in a 
quantum system that is still far from the classical limit. 

We have previously shown that quantum anomalous 
diffusion in a simple delta-kicked system can be faster 
than the underlying classical anomalous diffusion |30| . 



However, the spatial symmetry therein always gives a 
zero momentum expectation value and therefore one nec- 
essarily considers the momentum variance in order to 
make a meaningful quantum-classical comparison. By 
contrast, our ratchet accelerator model here allows for 
a study of the quantum-classical differences in the mo- 
mentum expectation value itself. As such, the role of 
the quantum tunneling between a chaotic sea and trans- 
porting regular islands becomes even more clear: while it 
can enhance the increase rate of the momentum variance 
, it suppresses the net acceleration rate of the average 
momentum in our ratchet accelerator model. 

It would be of great interest to consider a molecular re- 
alization of the ratchet accelerator, where the variable x 
is understood as an angle variable and p is understood as 
an angular momentum variable. Since a microwave field 
[l8Ll39| will create a kicking potential of the form cos(x) 
and an off-resonance laser field [Tfll EH ] can create a kick- 
ing potential of the form cos 2 (x) , a diatomic subject to 
both microwave and off-resonance laser fields, with their 
polarization direction perpendicular to one another, is a 
possible candidate for realizing the ratchet accelerator we 
propose. However, since in this case the diatomic should 
be at least described by a three-dimensional rigid rotor, 
which differs from the planetary kicked rotor in a num- 
ber of aspects [3^, detailed calculations should be car- 
ried out to guide the experimental studies. Nevertheless, 
we believe that this ratchet accelerator model may be of 
importance for manipulating the rotational motion of di- 
atomics and may provide a means of orienting molecules 
while accelerating rotation without using a biased exter- 
nal field. 
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Sciences and Engineering Research Council of Canada. 
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